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ABSTRACT 

We use the GALPROP model for cosmic-ray propagation to obtain a new estimate of the Galactic 
component of gamma rays, and show that away from the Galactic plane it gives an accurate prediction 
of the observed EGRET intensities in the energy range 30 MeV - 50 GeV. On this basis we re- 
evaluate the extragalactic gamma-ray background. We find that for some energies previous work 
underestimated the Galactic contribution at high latitudes and hence overestimated the background. 
Our new background spectrum shows a positive curvature similar to that expected for models of the 
extragalactic emission based on the blazar population. 

Subject headings: diffusion — cosmic rays — ISM: general — diffuse radiation — gamma rays: obser- 
vations — gamma rays: theory 



1. INTRODUCTION 

The extragalactic diffuse 7-ray background emission 
(EGRB) is a superposition of all unresolved sources 
of high energy 7-ray emission in the Universe. Ac- 
tive galactic nuclei (AGN) are the dominant class of 
7-ray emitters known to emit up to the highest ener- 
gies. There is a consensus that a population of un- 
resolved AGN may contribute to the EGRB, however 
predictions range from 25% up to 100%. Contribu- 
tions from other extragalactic sources have been p lau- 
sibly suggested: galaxy clusters l)Ensslin et alJUQQ^ . en- 
ergetic particles in the shock waves ass ociated with large 
scale cosmological s tructure formation lILoeb fc Waxmanl 
120001 iMiniatil 12002). or distant 7-ray burst events. Po- 
tentially, if reliably derived, the EGRB can also pro- 
vide important informat ion about the phase of baryon- 
antib aryon annihilation iJGao et al.lll990t IDolgpv fc Silk! 
1993), evaporation of primordial b lack holes (Hawking 
1974i iMaki. Mitsui, fc Oritolll99^1 . annihilation of so- 
called weakly interacting massive partic les (WIMPs) 
l).Iungman. Kamionkowski. fc Griestl 119 9^ and ext ra- 
galactic IR and optical photon spectra ( Steckerl ll999). 

The EGRB is the component of the diffuse emission 
which is most difficult to determine. Its spectrum de- 
pends much on the adopted model of the Galactic back- 
ground which itself is not yet firmly established. The 
isotropic, presumably extragalactic component of the 
diffuse 7-ray flux was first disc overed by the SAS-2 
satell i te and confirmed by EGR ET ({Thompson Sz Fichtell 
119821 iSreekumar et al.1 Il998j) . However, it is not 



correct to assume that the isotropic component is 
wholly extragalactic since even at the Galactic poles 
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it is comparable to the Galactic contribution from in- 
verse Compton s cattering of the Galactic plane pho- 
tons and CMB ([Strong. Moskalenko. & Reimerl 1200(1 
iMoskalenko fc Strong! I2000T1 . The determination of the 
EGRB is thus model-dependent and influenced by the 
adopted size of the Galactic halo, the electron spectrum 
there, and the spectrum of low-energy background pho- 
tons which must be derived inde pendently. 

Extensive work has been done (jSreekumar et alii 1998(1 
to derive the spectrum of the EGRB based on EGRET 
data. The relation of modelled Galactic diffuse emission 
to total measured diffuse emission was used to determine 
the EGRB, as the extrapolation to zero Galactic contri- 
bution. The derived spectral index —2.10 ±0.03 appears 
t o be close to t hat of 7-ray blazars. 

Dixon et al. (1998), using an model-independent ap- 
proach, also concluded that the EGRB will be affected 
by a significant contribution from a Galactic halo com- 
ponent. 

In a companion paper ijStrong. Moskalenko. &: Reimerl 
|2004|) we use the GALPROP code to infer a new model 
for Galactic diffuse continuum 7-rays. This model re- 
produces successfully diffuse 7-ray emission from the en- 
tire sky and gives a good linear prediction for observed 
vs. predicted 7-ray intensities. In view of the success of 
this model, we use it as the basis of a new determina- 
tion of the EGRB using EGRET data from 30 MeV to 
50 GeV. The GALPROP propa gation code and previous 
results are describ ed elsewhere IMoskale nko et af]l2002t 
iStrong et al.ll2004l and references therein). 

2. THE PROCEDURE 

2.1. EGRET data 

The details of the procedure of handling the EGRET 
data and convolution procedure are described in detail 
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elsewhere l|Strong et al J 1200^1 : here we provide a brief 
summary. 

We use the EGRET counts and exposure all-sky 
maps in Gal a ctic c oordinates with 0.5° binsize, as in 
iStrong et all (|200(1 . The sources of the 3EG cata- 
logue have b e en re moved by the procedure described in 
IStrong et alJ 1)2000(1 . fully consistent with the 3EG point 
source listings. Apart from the most intense sources, 
the removal of sources has little influence on the com- 
parison with models if sufficiently large sky segments 
are investigated. For the spectra, the statistical er- 
rors on the EGRET data points are very small since 
the regions chosen have large solid angle; the system- 
atic error dominates and we have adopted values in the 
range 10-30% depending on energy as sh o wn in Ta- 
ble H cf. 15% adopt ed by ISreekumar et alJ l|1998j) and 
lEsposito et alJ l|1999|) . The predicted model skymaps are 
convolved with the EGRET po int-spread function as de- 
scribed in IStrong et alJ l(2000() . Here we use additional 
EGRET data in the energy ranges 10-20, 20-50 and 50- 
120 GeV. Because the instrumental response of EGRET 
determined at energies above 10 GeV is less certain com- 
pared to energies below 10 GeV, it is required to account 
for additional systematic uncertainties. In particular the 
EGRET effective area can only be deduced by extrapo- 
lation from the calibrate d effective area at lower energies 
((Thompson et al.lll993al) . We accordingly adopt values 
of 0.9, 0.8, and 0.7 times the 4-10 GeV effective area, 
respectively. 

2.2. Optimized model for the Galactic diffuse emission 

In the companion paper ((Strong et al.ll2004() we com- 
pared a range of models of Galactic diffuse emission, 
based on our CR propagation code GALPROP, with data 
from the Compton Gamma Ray Observatory. There we 
exploit the fact that the models predict quite specific be- 
haviour for different sky regions and this provides a crit- 
ical test: the "correct" model should be consistent with 
the data in all directions. We show that a new model, 
with moderate changes of electron and nucleon spectra 
relative to the "conventional" model, can well reproduce 
the 7-ray data and is compatible with locally observed 
particle spectra considering the expected level of spatial 
fluctuations in the Galaxy. The 7-ray data comparisons 
were extended over the entire sky and to 100 GeV in 
energy. We also exploited the recent improved measure- 
ments of the local proton, Helium, as well as antiproton, 
and positron spectra which are used as constraints on the 
proton spectrum in distant regions. 

To fit the Galactic diffuse emissio n we used 7 test re - 
gions covering the sky (for details see lStrong et all 2004'l. 
The model uses the proton (and He) spectral shape at 
high energies derived from the local data. The secondary 
antiproton and positron data tracing the proton spec- 
trum on a large scale provide an important constraint on 
the intensity normalization of the average spectrum. The 
adopted average electron spectrum resembles the local 
one renormalized upwards by a factor of ^4 and consis- 
tent with synchrotron index measurements. The adopted 
electron and nucleon spectra are compatible with the di- 
rect measurement considering fluctuations due to energy 
losses and the stochastic character of cosmic ray sources 
and propagation. 

The optimized model fits the observed Galactic diffuse 



7-ray spectra in all test regions in the energy range 30 
MeV - 100 GeV. The proposed scenario implies a sub- 
stantial contribution from IC at all energies, but espe- 
cially below 100 MeV and above 1 GeV. Also IC dom- 
inates at latitudes \b\ > 10° at all energies. The agree- 
ment of longitude profiles with the EGRET data is gen- 
erally good considering that the model does not attempt 
to include details of Galactic structure. The agreement 
of latitude profiles with EGRET is also good, in partic- 
ular the reproduction of the high-latitude variation con- 
firms the importance of the IC component which is much 
broader than the gas-related 7r°-decay and bremsstrah- 
lung emission. The outer Galaxy latitude profiles are in 
excellent agreement with the data. 

2.3. Method 

Since despite its success the model is nevertheless not 
exact the best approach is to fit the observed intensi- 
ties with a free scaling factor, so that the EGRB is de- 
termined as the intercept, thus removing any residual 
uncertainty in the absolute level of the Galactic compo- 
nents. This is the same method as in ISreekumar et ail 
(1998), the difference lies in the model. To reduce the 
effects of Galactic structure, point sources etc. the 
fits are made excluding the plane, using the sky region 
360° < I < 0°, 10° < \b\ < 80°; ideally both IC and gas- 
related components would be left free but they are diffi- 
cult to separate statistically at high latitudes, so we make 
a linear fit to the total IC+7r°-decay-|-bremsstrahlung, 
with the scaling factor and EGRB as parameters. The 
fit and errors are based on a simple \ 2 analysis, with (/, b) 
bins 360° x 2° to obtain sufficient statistics (at least 10 
counts per bin were required). For comparison we also 
made fits to the entire sky; in this case IC and gas-related 
contributions are easily separated, so that fits with both 
components free can be made in addition to fits to the 
total. The fits in the two regions then give some indi- 
cation of the model-dependent systematic error in our 
EGRB estimates. 

3. DETERMINING THE EGRB 

The scaling factors determined for the fits in 360° < 
I < 0°,10° < |6| < 80° (Tabled) reflect the deviation s 
from the model, visible in Fig. 8 of IStrong et alJ {2004), 
and are typically between 0.85 and 1.15 which is satis- 
factory considering the EGRET systematic uncertainty 
~15% indicating a good agreement outside the plane. 
There is only one region (2-4 GeV) where the scaling fac- 
tor is as large as 1.35. For the all-sky fits the 1-2 and 2-4 
GeV factors are large in the case of the IC component, 
probably related to an underestimate of the ISRF or of 
the electron density. In the range 30-50 MeV, the gas- 
related factor is 1.31 which may indicate underestimated 
bremsstrahlung, but the EGRET data in this region are 
subje ct to a correction as a resu lt of instrumental calibra- 
tions ((Thompson et al]l!993bt s ee also discussion o f the 
Kniffen factor in Section 3.1 of IStrong et alJl2004() and 
thus more uncertain compared to other energies. These 
are the least reliable ranges of our EGRB determination 
reflecting th e discrepancy in the spectrum mentioned in 
Section 6 of IStrong etail <|2004fl . The EGRB is however 
not very sensitive to the scaling factor, since the inter- 
cept requires only a small extrapolation to zero Galactic 
flux. The observed versus (fitted) predicted intensities 
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Fig. 1.— Observed versus predicted intensities for region 360° < I < 0°,10° < \b\ < 80° (GALPROP model ID 500190). Predicted is 
Galactic only, intercept is EGRB. 



are shown in Fig. ^ (region \b\ > 10°) using the param- good accuracy and the EGRB well determined. 

eters in Tableland Table |2 These plots give an idea Table [21 presents the EGRB, integrated over each en- 

of the reliability of the analysis: the relation is linear to ergy bin since this is given directly by the fitting. The 
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Table 1. Scaling factors of model components corresponding to fits 
in Table El 



Energy, MeV 


10° < |6| < 80°, 
total 






All-sky 










Total 


1G 




Gas 




30 


-50 


i nnc -|_ n (VI 


1 


07 ±0.017 


\J •Z)iJ _1_ xJ.KJZj 


I 


31 ± 


16 


50 


70 


0.864 ± 0.02 





88 ±0.009 


0.66 ± 0.05 


1 


13 ±0 


06 


70 


100 


0.872 ±0.01 





84 ± 0.005 


0.82 ± 0.03 





86 ±0 


02 


100 


-150 


0.851 ±0.01 





82 ± 0.004 


0.85 ± 0.03 





80 ±0 


01 


150 


-300 


0.874 ±0.01 





83 ± 0.003 


0.95 ±0.02 





79 ±0 


01 


300 


-500 


0.965 ±0.01 





90 ± 0.004 


1.14 ±0.04 





84 ±0 


01 


500- 


1000 


1.02 ±0.01 


1 


09 ± 0.004 


1.28 ±0.04 





87 ±0 


01 


1000- 


-2000 


1.15 ±0.02 


1 


22 ± 0.007 


1.50 ±0.07 





98 ±0 


02 


2000- 


-4000 


1.35 ± 0.04 


1 


39 ±0.014 


1.54 ±0.11 


1 


11 ±0 


01 


4000- 


10000 


0.96 ± 0.07 





92 ± 0.022 


0.98 ±0.15 





89 ±0 


07 


10000 


-20000 


0.98 ± 0.04 





81 ±0.058 


0.62 ± 0.39 





91 ±0 


21 


20000- 


-50000 




0.87 ±0.16 


1.14 ± 1.3 





71 ±0 


73 



Note. — "Total" refers to fitting sum of all model diffuse components, "IC" and "gas" refer to fits separating IC from gas-related 
components. Only statistical errors are given. 

Table 2. Estimates of EGRB obtained by fitting optimized model 
500190 to EGRET data. 



10° < |5| < 80°, Adopted fractional All-sky 



Energy, MeV 


total 


systematic error 


Total 


lG±gas 


Sreekumar et al. ( 1998) a 


30-50 


16.8 ±0.66 


0.30 


16.0 ± 0.48 


17.0 ±0.84 


24.0 ±7.0 


50-70 


10.6 ±0.19 


0.15 


10.4 ±0.13 


11.1 ±0.22 


13.26 ±2.58 


70-100 


6.66 ±0.10 


0.10 


6.7 ± 0.065 


6.75 ±0.11 


7.83 ±1.05 


100-150 


4.48 ± 0.07 


0.10 


4.6 ± 0.045 


4.55 ± 0.076 


5.5 ±0.75 


150-300 


3.92 ± 0.06 


0.10 


4.2 ± 0.040 


3.92 ± 0.067 


5.4 ±0.72 


300-500 


1.20 ±0.04 


0.10 


1.39 ±0.025 


1.17 ±0.041 


1.97 ±0.268 


500-1000 


0.76 ± 0.04 


0.10 


0.93 ±0.023 


0.70 ± 0.037 


1.36 ±0.185 


1000-2000 


0.32 ± 0.03 


0.10 


0.39 ±0.018 


0.26 ±0.028 


0.617 ±0.084 


2000-4000 


0.24 ± 0.02 


0.10 


0.30 ±0.014 


0.25 ±0.022 


0.304 ± 0.044 


4000-10000 


0.117 ±0.02 


0.10 


0.13 ±0.010 


0.12 ±0.016 


0.1956 ±0.0288 


10000-20000 


0.025 ± 0.04 


0.25 


0.034 ±0.011 


0.04 ± 0.02 


0.053 ±0.016 


20000-50000 




0.25 


0.011 ±0.022 




0.0276 ± 0.0096 


50000-120000 










0.0147 ±0.0063 



Note. — "Total" refers to fitting sum of all model diffuse components, "IC" and "gas" refer to fits separating IC from gas-related 
components. For our fits statistical errors are given together with adopted systematic errors on the EGRET response. Units: 10 — 6 cm -2 
sr~ 1 s — 

a Values are from their Table 1 integrated over the energy bin. The error bars include systematic errors. 

values from lSreekumar et alJ l)1998f) (their Table 1) were 
thus integrated (assuming an E~ 2 spectrum) to compare 
with ours. It is seen that the two fitted regions (\b\ > 10° 
and all-sky) give consistent results, indicating that there 
is no large systematic effect; it shows a model-dependent 
systematic uncertainty of 5-25%. This is comparable to 
the ~15% systematic uncertainty on EGRET effective 
area. The adopted energy-dependent systematic error 
for EGRET data is also given in the Table. The total er- 
ror is obtained by combining the statistical and EGRET 
systematic errors in quadrature. 

Fig. [2] shows the spectrum of EGRB as derived in the 
present work (Table [3| . The data points of the differ- 
ential spectrum were obtained from Table [21 (columns 
1-3 and column 4 above 10 GeV) assuming a power-law 
spectrum within the energy bin. The two E > 10 GeV 
points were taken as given in "All-sky, Total" column, 
while the 20-50 GeV point was plotted as a la upper 
limit. Above 50 GeV the data were insufficient to give a 
significant result by our method. 

Our estimate s lie significantly below those of 
iSreekumar et all <|1998f) . in most energy ranges. The pos- 
itive curvature in our EGRB spectru m is interesting and 
not unexpected for a blazar origin ( Sala mon fc Steckerl 
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Fig. 2. — Comparison of our EGRB sp ectrum ( solid, red) as 
given in Table El with that from ISreekumar et alJ <1998T) (d ots, 
magenta). The dashed (blue) line is ttie TSreekuinar^t^ll fit: 
- 1 sr" 1 MeV" 1 . 
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Fig. 3. — Ex tragalactic X-ray and 7-ray spectrum. Data compilation from Srcckumar et al.l 1119981) except for COMPTEL 
IWeide nspointncr ct al. 2000) and EGRET 30 MeV - 20 GeV (this work). 



Table 3. EGRB intensity. 



Ej, MeV 


Intensity, 
cm" 2 sr _1 s _1 MeV -1 


Total error 


30-50 


8.40 x 10~ 7 


2.54 x 10~ 7 


50-70 


5.30 x 10- 7 


0.80 x 10~ 7 


70-100 


2.22 x 10" 7 


0.22 x 10~ 7 


100-150 


8.96 x 10" 8 


0.91 x 10" 8 


150-300 


2.61 x 10" 8 


0.26 x 10~ 8 


300-500 


6.00 x 10~ 9 


0.63 x 10- 9 


500-1000 


1.52 x 10" 9 


0.17 x 10" 9 


1000-2000 


3.20 x 10" 10 


0.44 x 10~ 10 


2000-4000 


1.20 x 10" 10 


0.16 x 10~ 10 


4000-10000 


1.95 x 10" 11 


0.30 x 10 -11 


10000-20000 


3.40 x 10" 12 


1.39 x 10~ 12 


20000-50000 


< 1.11 x 10~ 12 





Note. — EGRB intensity (as plotted in Fig. derived from 
columns 1—3 of Tableland column 4 above 10 GeV. 



199$ but in view of the systematic uncertainties should 
not be taken too liter ally; a similar, les s pronounced ef- 
fect is present in the iSreekumar et al.1 spectrum. Our 
spectrum is not consistent with a power law. The 0.1 
- 10 GeV intensity for the region outside the Galac- 
tic plane (360° < I < 0°, 10° < |6| < 80°) of 
(11.1 ± 0.1) x 10~ 6 cm~ 2 sr^ 1 s _1 can be compared with 
(14.5 ± 0.5) xlO -6 cm -2 sr _1 s" 1 from ISreekumar et al.1 



<199: 

Fig. [3] shows the extragalact ic X- and 7-ray back- 
ground, using the compilation by Sre ekumar et al.l l|1998f) 
but using our new EGRET values, and also u pdated 
COMPTEL results <)Weidenspointner et al.ll2000(l . 

4. FURTHER CHECKS FOR SYSTEMATICS DUE TO 
MODEL 

In order to check the robustness of our EGRB esti- 
mates, we compare estimates based on data in different 
hemispheres and four quarter spheres (Table 0J) . This 
tests for the presence of apparent anisotropies in the 
EGRB values, and thus gives an independent estimate 
of the systematic error due to model inadequacies. In 
order to have sufficient statistics we consider only the 
integral energy range 0.1-10 GeV. The hemisphere devi- 
ations are up to 16% from the fit to the full \b\ > 10° 
region. The quarter sphere values deviate by up to 25% 
from the full region. This may reflect essentially different 
exposures (e.g., southern hemisphere is less exposed than 
northern hemisphere). The deviations are also probably 
indicating lack of symmetry in the observed Galactic 7- 
ray emission not reflected in the model. However, using 
all regions together, as we have done for our results, will 
tend to compensate such systematics. 

Fig. 0| shows the EGRB determined from the fits to 
all the regions of Table 0] The quarter-sphere spectra 
show more scatter than the hemisphere spectra, as ex- 
pected since there are less data fitted. The extrema are 
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Table 4. Test for. model-dependent systematics. 



1 




b 


Intensity 0.1-10 GeV 


Description 


0-360 


< 


-10, > +10 


11.10 ± 0.12 


N+S hemispheres 


0-360 


< 


-10 


11.70 + 0.15 


N hemisphere 


0-360 


> 


+10 


9.28 + 0.21 


S hemisphere 


270-90 


< 


-10, > +10 


11.90 + 0.17 


Inner Galaxy N+S 


90-270 


< 


-10, > +10 


9.75 ± 0.17 


Outer Galaxy N+S 


0-180 


< 


-10, > +10 


10.80 + 0.17 


Positive longitudes N+S 


180-360 


< 


-10, > +10 


11.60 + 0.16 


Negative longitude N+S 


270-90 


> 


+10 


13.00 + 0.22 


Inner Galaxy N 


270-90 


< 


-10 


9.14 + 0.32 


Inner Galaxy S 


90-270 


> +10 


10.60 + 0.22 


Outer Galaxy N 


90-270 


< 


-10 


8.18 + 0.34 


Outer Galaxy S 



Note. — Units: 10 6 cm 2 sr 1 s 1 . Errors are statistical only. 




iff iff 
energy (MeV) 

Fig. 4. — EGRB spectrum determined for each of the regions 
in Table |3] Red bars: hemispheres, blue bars: quarter spheres. 
The upper and lower bound spectra are shown in green, and the 
adopted EGRB in black. E nergies are shifted slightly for clarity. 
The dashed (blue) line is the Srcckumar ct al. fit. 



plotted to give an indication of the extreme upper and 
lower bounds on the spectrum accounting for systematics 
due to the model. Compared with the adopted EGRB 
(also shown), the extreme upper limits lie 20-50% higher 
(100% for 4-10 GeV), the lower limits 10-50% lower. 



These extreme bounds are quite conservative; note that 
since this plot illustrates the model-dependent error, the 
additional error due to exposure uncertainty (Section 3) 
is not included here. The adopted spectrum using all the 
|6| > 10° data is expected to be the most robust against 
such systematics since large-scale asymmetry effects will 
tend to cancel. 



5. CONCLUSIONS 

Based on our new optimized model for the diffuse 
Galactic 7-ray emission, a new EGRB spectrum has 
been derived. I t is low er and steeper than found by 
iSreekumar et al.l l)1998j) : it is not consistent with a 
power-law, and shows some positive curvature as ex- 
pected for an origin in blazars. 
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